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ABSTRACT 


'The search for habitable environments and biomarkers in exoplanetary atmospheres is the holy grail 
of exoplanet science. The detection of atmospheric signatures of habitable Earth-like exoplanets is 
challenging due to their small planet-star size contrast and thin atmospheres with high mean molecular 
weight. Recently, a new class of habitable exoplanets, called Hycean worlds, has been proposed, 
defined as temperate ocean-covered worlds with Ho-rich atmospheres. Their large sizes and extended 
atmospheres, compared to rocky planets of the same mass, make Hycean worlds significantly more 
accessible to atmospheric spectroscopy with the JWST. Here we report a transmission spectrum of 
the candidate Hycean world, K2-18 b, observed with the JWST NIRISS and NIRSpec instruments in 
the 0.9-5.2 um range. The spectrum reveals strong detections of methane (CH4) and carbon dioxide 
(CO2) at 5e and 3c confidence, respectively, with high volume mixing ratios of ~1% each in a Hə- 
rich atmosphere. The abundant CH, and CO» along with the non-detection of ammonia (NH3) are 
consistent with chemical predictions for an ocean under a temperate H2-rich atmosphere on K2-18 b. 
The spectrum also suggests potential signs of dimethyl sulfide (DMS), which has been predicted to be 
an observable biomarker in Hycean worlds, motivating considerations of possible biological activity on 
the planet. The detection of CH, resolves the long-standing missing methane problem for temperate 
exoplanets and the degeneracy in the atmospheric composition of K2-18 b from previous observations. 
We discuss possible implications of the findings, open questions, and future observations to explore 


this new regime in the search for life elsewhere. 
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biology(74) — Biosignatures(2018) 


1. INTRODUCTION 


'The detection and characterisation of habitable-zone 
exoplanets is a major frontier in modern astronomy. Un- 
til recently, the quest for exoplanetary habitability and 
biosignatures have been focused primarily on rocky ex- 
oplanets, naturally motivated by the terrestrial expe- 
rience of life (Kasting et al. 1993; Meadows & Barnes 
2018). The extreme diversity of exoplanetary systems 
witnessed over the past three decades motivates consid- 
erations of new avenues in the search for life elsewhere. 
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Such an endeavour may open the doors to a wider range 
of habitable environments that may be more numerous 
and more favourable to atmospheric characterisation. 
Hycean worlds, a recently proposed class of habitable 
exoplanet, represent one such avenue that is accessible 
to current observational facilities (Madhusudhan et al. 
2021). 

Hycean worlds are a class of water-rich sub-Neptunes 
with planet-wide oceans underlying Hə-rich atmo- 
spheres. Such planets have a significantly wider hab- 
itable zone compared to terrestrial planets. With ex- 
pected radii between 1-2.6 Rg for masses between 1- 
10 Mg, Hycean planets represent a habitable subset of 
temperate sub-Neptunes that allow for a vast diversity 
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of atmospheric and internal structures (Madhusudhan 
et al. 2020, 2021; Piette & Madhusudhan 2020; Nixon 
& Madhusudhan 2021). Such planets are also poten- 
tially abundant in the exoplanet population given the 
predominance of exoplanets in the sub-Neptune regime 
(Fulton & Petigura 2018). The large volatile content in 
the interior of a Hycean world implies a lower density 
and, hence, larger radius and lower gravity, compared 
to a rocky planet of comparable mass. The low gravity 
and low atmospheric mean molecular weight (MMW) 
in turn result in a larger atmospheric scale height for 
a given temperature relative to terrestrial-like exoplan- 
ets with high-MMW atmospheres. These factors make 
Hycean worlds readily accessible for atmospheric charac- 
terisation, including potential biomarker detection, us- 
ing modest observing time with JWST (Madhusudhan 
et al. 2021; Phillips et al. 2021, 2022; Leung et al. 2022). 
The Hycean planet class was motivated by the demon- 
stration that the bulk properties of the habitable zone 
sub-Neptune K2-18 b (Montet et al. 2015; Cloutier et al. 
2017, 2019; Benneke et al. 2019a) are consistent with 
the possibility of a water-rich interior and a liquid water 
ocean at habitable temperatures and pressures under- 
lying a Ho-rich atmosphere (Madhusudhan et al. 2020). 
The planet has a mass of 8.63 + 1.35 Ma and radius 
of 2.61 + 0.09 Ra, with an equilibrium temperature of 
~250-300 K for an albedo between 0-0.3 (Cloutier et al. 
2019; Benneke et al. 2019a). While a Hycean interpreta- 
tion for K2-18 b is plausible and promising, a broad set 
of other internal structures and non-habitable surface 
conditions are also compatible with its bulk properties 
(Madhusudhan et al. 2020; Piette & Madhusudhan 2020; 
Nixon & Madhusudhan 2021), especially when consid- 
ering only cloud/haze-free atmospheres (e.g. Scheucher 
et al. 2020; Pierrehumbert 2023; Innes et al. 2023). 
Originally, the observed transmission spectrum of the 
planet in the near-infrared (1.1-1.7 um) with the HST 
WFC3 spectrograph suggested a Hə-rich atmosphere 
with strong H20 absorption (Benneke et al. 2019a; 
Tsiaras et al. 2019; Madhusudhan et al. 2020). How- 
ever, other studies highlighted the degeneracy between 
H2O and CH, in the observed HST spectrum (Blain 
et al. 2021; Bézard et al. 2022), and potential contribu- 
tions due to stellar heterogeneities (Barclay et al. 2021), 
rendering the previous H20 inference inconclusive. 
Atmospheric observations with JWST have the poten- 
tial to provide important insights into the atmospheric, 
surface, and interior conditions of K2-18 b. The planet 
has been theoretically demonstrated to be accessible 
to detailed atmospheric characterisation with a modest 
amount of JWST time, including the possibility of de- 
tecting prominent CNO molecules, such as H20, CH4, 


NHs, as well as several biomarkers, such as (CH3)2S 
or dimethyl sulfide (DMS), methyl chloride (CH3C1), 
OCS and others (Madhusudhan et al. 2021). The ma- 
jor molecules are expected to be detectable even in the 
presence of high-altitude clouds (Constantinou & Mad- 
husudhan 2022). Furthermore, several recent theoreti- 
cal studies have demonstrated that atmospheric abun- 
dances of prominent CNO molecules can be used to infer 
the presence of surfaces beneath H2-rich atmospheres in 
temperate sub-Neptunes (Yu et al. 2021; Hu et al. 2021; 
Tsai et al. 2021). For example, the presence of an ocean 
underneath a shallow H2-rich atmosphere, as would be 
the case for a Hycean world, may be inferred by an en- 
hanced abundance of CO2, H20 and/or CHy, but with 
a depletion of NH; (Hu et al. 2021; Tsai et al. 2021; 
Madhusudhan et al. 2023). 

In this work, we report the first JWST transmission 
spectrum of K2-18 b. The spectrum was observed us- 
ing NIRISS SOSS and NIRSpec G395H instruments in 
the 0.9-5.2 um wavelength range, which contains strong 
spectral features of multiple chemical species. The 
chemical constraints derived from the observed spec- 
trum provide key insights into its atmospheric and sur- 
face conditions and pave the way for a new era of at- 
mospheric characterisation of low-mass exoplanets with 
JWST. In what follows, we present our JWST observa- 
tions and data reduction in section 2. We discuss our 
atmospheric retrievals of the transmission spectrum in 
section 3. We summarize our results and discuss the 
implications in section 4. 


2. OBSERVATIONS AND DATA REDUCTION 


We report transmission spectroscopy of K2-18 b using 
the JWST NIRSpec (Ferruit et al. 2012; Birkmann et al. 
2014) and NIRISS (Doyon et al. 2012) instruments. We 
observed two primary transits of the planet in front of its 
host star, one with each instrument, as part of the JWST 
GO Program 2722 (PI: N. Madhusudhan). The first 
transit was observed using the NIRSpec G395H grating 
between Jan 20, 2023, 18:37:38 UTC and Jan 21, 2023, 
01:11:32 UTC for a total exposure time of 5.3 hours, 
which is nearly twice the expected transit duration. The 
observation was made in the Bright Object Time Se- 
ries (BOTS) mode with the F290LP filter, the SUB2048 
subarray and the NRSRAPID readout pattern, with the 
spectra dispersed over two detectors (NRS1 and NRS2). 
The two detectors, NRS1 and NRS2, span wavelength 
ranges of 2.73-3.72 um and 3.82-5.17 um, respectively, 
with a gap in between at 3.72-3.82 um. The G395H 
grating offers the highest resolution mode of NIRSpec 
with R~2700. The spectroscopic time-series observa- 
tion comprised of 1625 integrations, with 12 groups per 
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integration. For NIRSpec, the host star K2-18 was too 
bright for target acquisition (TA). Therefore, another 
nearby target (2MASSJ11301306+0735116) within the 
splitting distance of the science target was used for TA. 

The second observation was conducted using the 
NIRISS Single Object Slitless Spectroscopy (SOSS) in- 
strument mode (Albert et al. 2023) between Jun 1, 2023, 
13:49:20 UTC and Jun 1, 2023, 19:36:05 UTC, totalling 
an exposure time of 4.9 hours. The observation used the 
GR700XD grism (R~700), the CLEAR filter, the SUB- 
STRIP256 subarray and the NISRAPID readout pat- 
tern, giving a wavelength coverage of 0.85- 2.85 um for 
the first spectral order. The exposure consisted of 648 
integrations, with 4 groups per integration. There were 
no tilt events or high-gain antenna movements during 
any of the observations. 


2.1. NIRSpec 


'The data reduction is conducted using a combination 
of the JWST Science Calibration Pipeline (Bushouse 
2020) and our custom-built pipeline for the spectral ex- 
traction. We start with the raw 2D images (the .uncal 
files) of the spectroscopic time-series data which contain 
group-level counts for each integration. Stage 1 of the 
data reduction is performed mainly using the JWST Sci- 
ence Calibration Pipeline. This involves performing sat- 
uration flagging, superbias subtraction, reference pixel 
correction, linearity correction, dark current subtrac- 
tion (where reference data is available), jump detection 
and linear fitting of the group-level ramps to obtain the 
count rate for each pixel per integration. This is re- 
peated for each integration in the exposure. The jump 
detection threshold is set at 5o. Prior to ramp fitting, we 
also perform an additional step for background subtrac- 
tion, at the group level, in order to mitigate 1/f noise, 
as common in previous works (JWST Transiting Exo- 
planet Community Early Release Science Team et al. 
2023; Rustamkulov et al. 2023; Alderson et al. 2023). 
For this step we measure the background level for each 
detector column, as the mean of the pixels that are +10 
pixels away from the mid-point of the curved trace, while 
masking bad pixels and cosmic ray hits. The outputs of 
Stage 1 are 2D images of the count rate for each inte- 
gration, saved as .rateints files. 

We also use steps from the JWST Science Calibration 
Pipeline for Stage 2, which applies the wavelength cali- 
bration for the spectral trace. Following previous studies 
using NIRSpec (e.g. Alderson et al. 2023), we forego the 
flat field correction in Stage 2 for our differential tran- 
sit measurement. The resultant 2D images along with 
the wavelength calibration, which are saved as .calints 


files, are then used for the spectral extraction of the time 
series of 1D stellar spectra. 

We conduct the spectral extraction, Stage 3, applying 
our custom-built pipeline to convert the 2D images into 
1D spectra. This is conducted for each of the two detec- 
tors (NRS1 and NRS2) separately. We first create a bad 
pixel mask based on the data quality flags in the Stage 
2 products. We then extract the 1D spectrum from the 
2D image using an optimal extraction algorithm (Horne 
1986). To perform the extraction, we obtain the point- 
spread functions (PSFs) as the sum of the first three 
principal components of the time series of the detector 
images!, inspired by the principal component analysis 
based morphology analysis in Coulombe et al. (2023). 
This takes into account the wavelength and time de- 
pendence of the PSF. Outliers were iteratively rejected 
during the spectrum extraction, with the threshold set 
at 5c. Spectral channels with more than 20% of the flux 
masked were discarded from further analysis. 


2.2. NIRISS 


'To reduce the NIRISS data, we used the JWST Sci- 
ence Calibration Pipeline (Bushouse 2020) for Stages 1- 
2 and the JExoRES pipeline for the spectral extraction 
(Holmberg & Madhusudhan 2023). During Stage 1, we 
perform the standard saturation flagging, superbias sub- 
traction, linearity correction, jump detection (threshold 
set at 5c) and fitting of the group-level ramp to obtain 
the count rate. We perform a custom background sub- 
traction step before the linearity correction and ramp 
fitting to reduce the effect of 1/f noise. In line with 
Radica et al. (2023); Albert et al. (2023), this involves 
temporarily subtracting a model of the flux of the de- 
tector, containing both the background and stellar flux, 
to reveal the 1/f noise. Initially, we use the group-wise 
median stack of all integrations to model the flux. For 
each group and integration, the 1/f noise level is then 
estimated using the median of each column, while mask- 
ing the traces of the spectral orders as well as bad pixels. 
This 1/f noise estimate is subtracted from the data to 
reduce the level of correlated noise. We later repeat 
all stages of the NIRISS data reduction using updated 
models of the background level and stellar flux to refine 
the 1/f noise correction step, as described further be- 
low. Furthermore, we do not perform the dark current 
subtraction, given the quality of the dark reference file 
from the Calibration Reference Data System (CRDS) 
is insufficient, so as not to contribute additional noise 
(Feinstein et al. 2023; Radica et al. 2023). 


1 For this the sum of the detector columns was normalised to unity 


and bad pixels were interpolated. 
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Figure 1. Left: White light curve for the transit of K2-18 b observed with NIRSpec G395H. The top panel shows the combined 
white light curve from NRS1 and NRS2, together with the lo model interval in red. The bottom panel shows the residuals 
after subtracting the median model, together with line at zero. The standard deviation of the residuals is measured to be 123 
ppm, which is 1.2x the expected noise level. The red line is shown to indicate zero. Right: The top panel shows the normalised 
spectroscopic light curves (binned in wavelength for visual clarity). The detector gap is shown in white. The bottom panel 
depicts the transit depth precision at pixel resolution. 
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Figure 2. Left: White light curve for the transit of K2-18 b observed with NIRISS SOSS. The blue contour shows the lo 
model interval. We find evidence of a spot occultation at the start of the transit. The bottom panel shows the residuals after 
subtracting the median model. The standard deviation of the residuals is measured to be 115 ppm, which is 2.2x the expected 
noise level. The blue line is shown to indicate zero. Right: The top panel shows the normalised spectroscopic light curves 
(binned in wavelength for visual clarity). The bottom panel shows the transit depth precision obtained from the light curve 


fitting, binned at two pixel columns per light curve. 


For Stage 2, we perform the flat field correction be- 
fore modelling the background flux in a two-step process. 
First, we use the background model available via JDOX, 
created from program 1541, scaled to the median of all 
integrations in a small rectangular region (x € [720, 770], 
y € [210,250]). By subtracting this we precisely remove 
the brighter of the two background components, caused 
by zodiacal light. However, this leaves residuals in the 
dimmer component of the background for columns up to 
the —700 pixel column. This is corrected by first gen- 
erating a median image from all integrations, and then 
from this median image obtaining the median of each 
pixel column (while masking spectral traces, contami- 
nation and bad pixels). These column-wise medians are 
then subtracted from each integration image. This is 


only performed for columns up to the 700 pixel column. 
We use this additional background flux to update the 
(scaled) background model. Next, we perform the order 
tracing, PSF estimation, and spectral extraction (Stage 
3) according to Holmberg & Madhusudhan (2023), with 
an extraction aperture of 35 pixels, leaving out the back- 
ground refinement since we correct the 1/f noise at the 
group level. As in the case of NIRSpec, spectral channels 
with more than 2096 of the flux masked were discarded 
from further analysis. 

Finally, we repeat the data reduction stages in order 
to improve the 1/f noise correction. This time, we model 
the flux at the group level using the updated background 
model from above and the group-wise out-of-transit me- 
dian stellar flux. The background model is scaled us- 
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ing the integration time of each group while the out- 
of-transit stellar flux (minus the background) is scaled 
with a model of the transit light curve, derived using 
the initial white light curve from the first spectral order. 
We note that in the remaining analysis we only consider 
the first spectral order given that the second order has 
a considerably lower flux level, meaning that it is more 
sensitive to systematics due to sources of contamination. 


2.3. Starspot occultations 


Starspot and faculae crossings during exoplanet tran- 
sits are known to affect the apparent transit depth (e.g. 
Pont et al. 2008; Czesla et al. 2009). If left uncorrected, 
starspot and faculae occultations effectively decrease 
and increase the observed transit depth, respectively. 
'This effect is wavelength-dependent and must therefore 
be accounted for so as not to impact the transmission 
spectrum. From the present observations we find evi- 
dence of starspot occultations, which is especially strong 
for the transit observed with NIRISS. Fortunately, cor- 
recting the transmission spectrum is possible by measur- 
ing the intensity ratio and size of the occulted feature 
from the transit light curves. 

We perform a joint inference of the transit parameters 
and the properties of the active region using the semi- 
analytical spot modelling code SPOTROD (Béky et al. 
2014). Hence, we also avoid having to discard the af- 
fected data. The code computes the transit light curve 
with arbitrary limb darkening (affecting the stellar pho- 
tosphere and spots equally) and homogeneous circular 
star spots or faculae. The spot is represented by four pa- 
rameters: the spot-star radius ratio Rspor/R., the spot 
to unspotted stellar surface intensity ratio f, and the 
coordinates of the spot centre projected onto the stellar 
surface (0, r?). An intensity ratio below or above one 
represent starspots or faculae, respectively. Similar to 
Espinoza et al. (2019), we use nested sampling (Skilling 
2004), implemented via MultiNest (Feroz et al. 2009), 
to obtain the Bayesian model evidence and parameter 
estimation. This allows us to perform a Bayesian model 
comparison between transit models with or without stel- 
lar spots. For this, we use uniform priors between 0 and 
0.3 for the spot-star radius ratio, 0 and 2 for the in- 
tensity ratio, 0 and 27 for ? (with periodic boundary 
condition), and 0 and 1 for r?, ensuring uniform sam- 
pling of the disk. 

Using the NIRISS white light curve from the first or- 
der, we find that a single starspot is strongly preferred 
over no spots with a Bayes factor of In B — 21.6, cor- 
responding to a significance of 6.90. Fig. 2 shows the 
white light curve and the fitted model. We obtain an 
intensity ratio of f = 0.9329*0005; and a spot-star ra- 


dius ratio of Rspot/Rs = 0.254* 0041. For the NIRSpec 
white light curve (NRS1 and NRS2 combined), includ- 
ing a spot is also marginally preferred by the data com- 
pared to a model without spots, with a Bayes factor 
of In B = 1.15, corresponding to a significance of 2.1c. 
'The spot parameters in this case are f — 0350. and 
Rspot/ Be = 0.115100. For these reasons, we choose 


to use the spot modelling for both observations. 


2.4. Light curve analysis 


The 1D spectral time series from both observations 
are then used for light curve fitting to derive the transit 
depths. This is done in three stages. The first stage 
uses the white light curves, from both observations, to 
derive the wavelength-independent system parameters 
at high precision. For NIRISS we use the white light 
curve from the first order and for NIRSpec we use the 
combined white light curve from both NRS1 and NRS2. 
In the next stage, we bin the light curves to R ~ 20 
and fit the wavelength-dependent limb-darkening coef- 
ficients (LDCs). Finally, we fix these LDCs in the re- 
spective R ~ 20 bins and fit the transit depths at native 
resolution to obtain the final transmission spectrum of 
the planet. 

As described above, we model the transit light curves 
using SPOTROD (Béky et al. 2014) as we find evidence 
of starspot occultations. We assume a circular orbit 
with the orbital period from Benneke et al. (2019a). We 
adopt the two-parameter quadratic limb-darkening law, 
in line with previous JWST transmission spectroscopy 
studies of M dwarf systems (e.g. Moran et al. 2023; 
Lustig-Yaeger et al. 2023) and previous work on K2- 
18 b (Benneke et al. 2019a). For the LDCs, we use 
the parameterization and priors by Kipping (2013). To 
model the baseline flux we implement both a linear and a 
quadratic trend. We find that the NIRISS out-of-transit 
white light curve shows a preference for a linear trend, 
while the NIRSpec data show a weak preference for a 
quadratic trend in the case of NRS1 and a moderate 
preference in the case of NRS2. For the white light 
curves, we use a linear trend for both instruments, given 
that we find the derived system parameters to be insen- 
sitive to the choice of trend. However, for the spectrum, 
we consider both trends for NIRSpec. For NIRISS, we 
fix the trend to be linear. Moreover, we discard the 
first 5 minutes of both observations due to a small set- 
tling ramp. Apart from starspot occultations, no other 
systematics can be identified, highlighting the excellent 
data quality provided by JWST. 

For the white light curve of each observation (see left 
panels of Figs. 1 and 2), we fit for To, the mid-transit 
time, a/ R., the normalised semi-major axis, i, the or- 
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Figure 3. The transmission spectrum of K2-18 b. The observed JWST spectrum and retrieved model fits are shown for the 
l-offset retrieval case discussed in section 3. The data in orange show our NIRISS spectrum between 0.9 - 2.8 um and those 
in dark red show our NIRSpec G395H spectrum between 2.8 - 5.2 um. The spectra are binned to R = 25 and R = 55 for 
NIRISS and NIRSpec, respectively, for visual clarity. The retrievals are conducted on the native resolution spectra, resulting in 
the best-fit reduced x2 = 1.080. The NIRSpec spectrum is vertically offset by —41 ppm, corresponding to the median retrieved 
offset in the 1-offset retrieval case. The blue curve shows the median retrieved model spectrum, while the medium and lighter 
blue contours denote the lo and 2c intervals, respectively. Yellow points correspond to the median spectrum binned to match 
the observations. The prominent molecules responsible for the features in different spectral regions are labelled. 


NIRISS 
60096.7293687) 000062 


Parameter 
Mid-transit time, To (BJD - 2400000.5) 


NIRSpec 
59964.969453. 9.000085 


Weighted average 


Inclination, i (°) 89.5501 0-021 89.567+9:012 89.563 + 0.010 
Normalised semi-major axis, a/R. 79.9*14 80.92 0-75 80.68 + 0.68 
Planet-to-star radius ratio, Rp/R. 0.05412+9:90019 0.05441Ż0:00021 - 

First LDC, ui 0.06540: EE 0.301 Fire - 


Second LDC, u» 0.220+9:966 —0.098* 0-097 
Table 1. Parameter estimation from the white light curve analysis of our JWST NIRISS SOSS and NIRSpec G395H observations 
of K2-18 b in the 0.9-5.2 um wavelength range. For the white light curve of NIRSpec, we combine NRS1 and NR$S2. For both 


white light curves, we use a linear trend to model the baseline flux as well as the spot modelling described in section 2.3. The 
orbital period is held fixed at P = 32.940045 days, adopted from Benneke et al. (20192). 


bital inclination, R,/R,, the planet-to-star radius ratio, 
uj; and uz, the quadratic LDCs, two parameters rep- 
resenting the baseline flux, and four spot parameters 
(described in 2.3). In the likelihood, we also include 
a parameter to inflate the photometric uncertainties to 
match the residual scatter between the data and the 
transit light curve model. We measure the precision 
to be 1.2x and 2.2x the expected noise level (photon 
and read noise, propagated using the Jacobian of the 
transit model) for the white light curves from NIRSpec 
(NRS1 and NRS2 combined) and NIRISS, respectively. 
To sample the posteriors and estimate the parameters 
we use MultiNest (Feroz et al. 2009), resulting in the 
parameters given in Table 1. For all cases, we find 
that the derived system parameters are consistent within 


1c, regardless of the trend and the choice to model the 
starspot or not. 

Next, we bin the spectroscopic light curves to R ~ 
20 for the purpose of fitting the wavelength-dependent 
LDCs. This resolution strikes a balance between preci- 
sion and the expected wavelength variability of the limb 
darkening. We choose to fit the LDCs instead of us- 
ing values from stellar atmospheric models in order to 
maximise accuracy (Csizmadia et al. 2013; Espinoza & 
Jordán 2015). We fit these binned light curves using 
MultiNest and fix the system parameters (To, a/Rx, i) 
to the values obtained from the white light curve analy- 
sis (see Table 1), i.e. the weighted average a/R,, i, and 
To from each observation. We also fixed Rspot/ R», 0, 
and r? to the best-fit parameters for each white light 
curve. Equipped with empirical LDCs, we go on to fit 
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the transit depths of the high-resolution light curves, 
while fixing the LDCs to the values within their re- 
spective R ~ 20 bin. This leaves only Ry/R., f, the 
2-3 trend parameters, and the uncertainty scaling pa- 
rameter as free parameters. We fit the light curves at 
the pixel level for NIRSpec and two pixels per bin for 
NIRISS, given the potential inaccuracy of the NIRISS 
SOSS wavelength calibration (Albert et al. 2023). To 
fit these high-resolution light curves we use the Leven- 
berg—Marquardt algorithm, as applied in previous works 
(Alderson et al. 2023; Moran et al. 2023). The right 
panel of Figs. 1 and 2 show the precision of the result- 
ing transmission spectrum at high resolution. In total, 
we have 1010 and 3401 spectral data points for NIRISS 
and NIRSpec, respectively, covering the 0.9 - 5.2 um 
wavelength range. 

As mentioned above, we obtain two spectra in the case 
of NIRSpec, one with a linear trend and another with 
a quadratic trend. For the quadratic trend case, we 
construct white light curves for NRS1 and NRS2, and 
fit these separately to obtain detector-averaged values 
for the quadratic trend component?. For each detector, 
we then fix these values when fitting the spectroscopic 
light curves (Moran et al. 2023). Overall, we find that 
the NRS1 spectrum is almost agnostic to the choice of 
trend (~10 ppm difference), whereas, for NRS2, there 
is a significant difference between the two spectra (ap- 
proximately a 60 ppm offset). Since we only have one 
transit observation for NIRSpec, we use an offset as a 
free parameter in the atmospheric retrieval to account 
for potential baseline shifts. 


3. ATMOSPHERIC RETRIEVAL 


The observed transmission spectrum allows us to re- 
trieve the atmospheric properties of K2-18 b at the day- 
night terminator region. We perform the retrieval us- 
ing the AURA retrieval code (Pinhas et al. 2018) fol- 
lowing a similar approach to previous retrieval studies 
of the planet considering HST and/or simulated JWST 
observations (e.g. Madhusudhan et al. 2020, 2021; Wel- 
banks et al. 2019; Constantinou & Madhusudhan 2022). 
The planet’s terminator is modelled as a plane-parallel 
atmosphere in hydrostatic equilibrium, with uniform 
chemical composition. The chemical abundances and 
pressure-temperature (P-T) profile are free parameters 
in the model. The retrieval framework follows a free 
chemistry approach, whereby the individual mixing ra- 


2 We model the baseline flux as F(t) = Foyt(1+pit+ pat”), where 
Fout is the out-of-transit flux at the start of the observation, 
pi is the linear trend parameter, and p2 is the quadratic trend 
parameter. 


tio of each chemical species is a free parameter. The 
atmospheric temperature structure is modelled with a 
parametric P-T profile (Madhusudhan & Seager 2009) 
with six free parameters. The model also considers in- 
homogeneous clouds/hazes at the day-night terminator 
region (MacDonald & Madhusudhan 2017; Pinhas et al. 
2018). They are modelled as a combination of a grey 
cloud deck at a parametric cloud top pressure P., above 
which are hazes with Rayleigh-like spectral contribu- 
tions, with an enhancement factor a and a scattering 
slope y; for Hə Rayleigh scattering, a = 1 and y = -4. 
The combined clouds and hazes cover a fraction ¢ of the 
atmosphere at the terminator region. 

The model includes molecular opacity contributions 
from prominent CNO molecules expected in temper- 
ate H5-rich atmospheres as considered in previous works 
(Pinhas et al. 2018; Welbanks et al. 2019; Constantinou 
& Madhusudhan 2022). These include H20, CH4, NH3, 
HCN, CO and CO2. The molecular absorption cross- 
sections are obtained following recent works (Gandhi & 
Madhusudhan 2017; Gandhi et al. 2020) using line lists 
from the following sources: H2O (Polyansky et al. 2018), 
CH, (Hargreaves et al. 2020), NH3 (Coles et al. 2019), 
HON (Harris et al. 2006; Barber et al. 2014), CO (Roth- 
man et al. 2010; Li et al. 2015) and CO» (Huang et al. 
2013, 2017). Pressure broadening due to Hə is consid- 
ered for all these molecules as described in Gandhi et al. 
(2020). 

We additionally consider five molecules which have 
been suggested to be promising biomarkers in habitable 
rocky exoplanets (Segura et al. 2005; Domagal-Goldman 
et al. 2011; Seager et al. 2013a,b; Catling et al. 2018; 
Schwieterman et al. 2018) as well as Hycean worlds 
(Madhusudhan et al. 2021): (CH3)3S or dimethyl sul- 
fide (DMS), CS2, CH3Cl, OCS and N20. The absorp- 
tion cross-sections of CH3Cl, OCS and N2O were com- 
puted from the corresponding line lists from the HI- 
TRAN database (Gordon et al. 2017): CH3Cl (Bray 
et al. 2011; Nikitin et al. 2016), OCS (Bouanich et al. 
1986; Golebiowski et al. 2014; Miiller et al. 2005; Auwera 
& Fayt 2006; Sung et al. 2009; Toth et al. 2010; Régalia- 
Jarlot et al. 2002), and N3O (Daumont et al. 2001). As 
pressure broadening due to H» is not available for these 
molecules, we only use thermal broadening. For DMS 
and CSz2, we use the absorption cross-sections provided 
directly by HITRAN (Kochanov et al. 2019; Sharpe 
et al. 2004; Gordon et al. 2017) at 1 bar and 298 K, fol- 
lowing Madhusudhan et al. (2021). In addition to molec- 
ular cross-sections, we also consider H2-H2 and H2-He 
collision-induced absorption (Borysow et al. 1988; Orton 
et al. 2007; Abel et al. 2011; Richard et al. 2012). 
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Figure 4. Retrieved posterior probability distributions for the mixing ratios of important molecules for the three retrieval 
cases described in section 3. The horizontal error-bar in the top row denotes each distribution’s median and corresponding lo 
interval for the three molecules with significant detections. The arrows in the lower row indicate 95% upper-limits. We find 
strong evidence for CH4 and CO», at a significance of 50 and 30, respectively. We find marginal evidence for DMS, and no 
significant evidence for the remaining molecules. The abundance estimates and detection significances are shown in Table 2. 


Cases CHa CO2 DMS H20 NH; CHCl CO HCN 
No offset —2.04*051 (4.70) —1.75t953 (2.90) —446*01 (2.40) «—3.21 «—446 <-—2.50 <-3.00 <-—2.41 
1 offset | —1.74*0 55 (5.00) —2.0916:31 (2.90) —6.35*158 (~1lo) <-3.06 «—4.51 <-3.80 «—3.50 «—2.92 


2 offsets  —1.89*055 (5.00) —2.05*020 (3.20) —6.87!15: (—-) | «—3.49 <-493 <-3.62 «-319 <-3.21 


'Table 2. Retrieved molecular abundances and detection significances of prominent molecules in the atmosphere of K2-18 b. The 
three canonical atmospheric model cases are described in section 3, and pertain to different considerations for offsets between 
data from different instruments. The molecular abundances are shown as log; of volume mixing ratios. The retrieved median 
and lo estimates are given for CH4, CO2 and DMS which show strong to marginal detections, and 95% upper-limits are given 
for the remaining molecules. The quantities in brackets for CHa, CO2 and DMS show the detection significances greater than 
lo. The detection significances have a nominal statistical uncertainty of ~0.1ø due to the uncertainty on the Bayesian evidence 
estimated by the nested sampling algorithm. See section 3.4. 


3.1. Retrieval Setup pling algorithm (Feroz et al. 2009) implemented through 
PyMultiNest (Buchner et al. 2014). The retrieval set- 
up and priors on the model parameters are similar to 
those in recent implementations of the AURA retrieval 
framework (Madhusudhan et al. 2020; Constantinou & 
Madhusudhan 2022) and are shown in Appendix C. We 
also consider variations to the canonical model, includ- 
ing a cloud/haze-free atmosphere, Mie scattering due to 
hazes, and the presence of stellar heterogeneities influ- 


Our canonical model comprises of 22 free parameters 
overall: 11 corresponding to the individual mixing ra- 
tios of the above chemical species, 6 for the P-T profile, 
4 for the clouds/hazes and 1 for the reference pressure 
Per, defined as the pressure at a fixed planetary radius 
of 2.61 Ra. The Bayesian inference and parameter es- 
timation is conducted using the MultiNest nested sam- 
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encing the spectrum, as discussed below. We use the 
present JWST NIRISS and NIRSpec transmission spec- 
tra of K2-18 b in the 0.9 - 2.8 um and 2.8 - 5.2 um 
ranges, respectively, at their native resolution for the 
retrieval. As discussed in section 2.4, we derive two spec- 
tra for NIRSpec corresponding to the two trends (linear 
and quadratic) we use to model the baseline flux; the 
NIRISS spectrum has a robust preference for a linear 
trend. Furthermore, the NIRSpec G395H grating uses 
two detectors (NRS1 and NRS2) and recent studies have 
suggested the possibility of an offset in the spectrum de- 
rived from a given detector (Moran et al. 2023). There- 
fore, we conducted retrievals on different combinations 
of NIRSpec spectra obtained with different trends and 
offsets. 

We consider two broad combinations of data: (1) 
NIRISS and linear-trend NIRSpec spectra and (2) 
NIRISS and quadratic-trend NIRSpec spectra. For each 
combination, we consider a range of different spectral 
offsets as free parameters in the retrieval. We consider 
four cases: a baseline case with no offsets, one com- 
bined offset for the NIRSpec spectrum, one offset for 
the NIRISS spectrum, and two separate offsets for the 
NIRSpec NRS1 and NRS2 spectra. The last case is the 
most conservative and is motivated by the transit depth 
offset between NRS1 and NRS2 recently reported by 
Moran et al. (2023). The offset on either NIRISS or 
NIRSpec represents cases where we assume no offset be- 
tween NIRSpec NRS1 and NRS2 (Alderson et al. 2023; 
Lustig-Yaeger et al. 2023; August et al. 2023) but in- 
stead allow for an offset between NIRISS and NIRSpec 
as a whole. In the no-offset case, we consider the data 
as is and do not perform any offsets. We conduct the 
four retrievals for each data combination and assess their 
relative Bayesian evidence. 

By comparing the Bayesian evidence of the considered 
cases, we find that generally some offset is preferred over 
no offsets. For the linear trend data, a single offset be- 
tween the data sets (either on NIRISS or NIRSpec) is 
the most preferred. The two one-offset cases are com- 
parable, with an offset on NIRSpec being marginally 
favoured over that on NIRISS. For the quadratic trend 
data, the two offset case is the most favoured. The fact 
that the configuration with separate offsets for NRS1 
and NRS2 is strongly preferred for the quadratic trend 
while not being favoured for the linear trend suggests 
that, for the present observation, using a quadratic trend 
contributes to an offset between NRS1 and NRS2. Based 
on these findings, from across the different combinations 
we finally select three nominal cases for NIRSpec along 
with NIRISS: (a) NIRSpec with a linear trend and no 
offsets, representing the data without modification, (b) 


NIRSpec with a linear trend and one offset, and (c) NIR- 
Spec with a quadratic trend and two separate offsets 
for NRS1 and NRS2, representing the most conserva- 
tive case. These are the retrieval cases considered in the 
rest of this work. In what follows we report the atmo- 
spheric properties at the day-night terminator region of 
K2-18 b retrieved using the transmission spectrum in 
each of these cases, as shown in Tables 2 and 3. 


3.2. Prominent CNO Molecules 


Our atmospheric retrieval provides important con- 
straints on the dominant CNO molecules expected in 
H2-rich atmospheres. The retrieved spectral fit is shown 
in Fig. 3, and the corresponding posterior distributions 
for several molecules are shown in Fig. 4. Amongst 
the prominent CNO molecules, we find strong spec- 
tral contributions from CH4 and CO» in a Hs-rich at- 
mosphere. For our retrieval with no offset, we derive 
log volume mixing ratios of log(Xcg,) = —2.04*09) 
and log(Xco,) = —1.75*0 03. For the one-offset case, 
we obtain log(Xcu,) = —1.74*0 59 and log(Xco,) = 
—2.09*031. As shown in Table 2, these abundance esti- 
mates for both molecules are consistent across the three 
retrieval cases, with median abundances of ~1% and 
average uncertainties below 1-dex, underscoring the ro- 
bustness of the derived estimates. Both CH4 and CO% 
are detected for the first time in a sub-Neptune exo- 
planet and the precision of their abundance estimates 
are the best measured for any molecule in a sub-Neptune 
atmosphere to date. 

We do not find significant contributions due to H20 
or NH3, but find 95% upper limits of -3.21 for log(Xu,o) 
and -4.46 for log(Xnu,) in the no-offset case. These up- 
per limits are also consistent with those from the other 
retrieval cases, as shown in Table 2. The non-detections 
of both molecules are important considering their strong 
spectral features and detectability expected in the 0.9- 
5.2 um range (Madhusudhan et al. 2021; Constantinou 
& Madhusudhan 2022). The non-detection of H20 is at 
odds with its previous inference using the HST WFC3 
spectrum in the 1.1-1.7 uim range (Tsiaras et al. 2019; 
Benneke et al. 2019a; Madhusudhan et al. 2020). A 
strong degeneracy between H20 and CH, in the HST 
WFC3 band was noted previously (Blain et al. 2021; 
Bézard et al. 2022). Our retrieved CH4 abundance is 
consistent with previous predictions of stronger absorp- 
tion due to CH, relative to H20 in the HST WFC3 
band (Blain et al. 2021; Bézard et al. 2022) and some 
upper bounds on the CH4 abundance (Madhusudhan 
et al. 2021; Blain et al. 2021; Bézard et al. 2022). Our 
results, therefore, resolve the degeneracy in the atmo- 
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Cases Tiombar (K) @ log a y log(P./bar) log(P.et/bar) OS1/ppm OS2/ppm x2 

No offset 257777" 0.63022 7.31433% «11.6773 95  —0.55*098 | —1.88*0$5 - : 1.082 
1 offset 242*72 0.087012 8.207123 -—11.11*291  —0.51'0998 2.4818 47 -41*H - 1.080 
2 offsets 41575 0.644913 8.214125 11.344233 —0.46t995 A007. — 3075 2619 1.081 


Table 3. Retrieved temperature, cloud/haze properties and reference pressure for the atmosphere of K2-18 b, as well as 
retrieved dataset offsets where applicable. Similarly to Table 2, the three canonical atmospheric model cases are described in 
section 3, which pertain to different considerations for offsets (denoted OS above) between data from different instruments. 
The 1-offset retrieval considers shifts to the NIRSpec observations relative to NIRspec, while the 2-offset retrieval considers 
two shifts, applied to observations from the NIRSpec NRS1 and NRS2 detectors. The temperature constraints shown are at 
10 mbar, which corresponds to the observed photosphere. In all cases, the best-fit reduced x2 is close to unity, with the degrees 
of freedom being 4389, 4388, and 4387 for no, 1, and 2 offsets, respectively. 


spheric composition of K2-18 b from previous observa- 
tions. 

While our NIRISS spectrum is generally in good agree- 
ment with the previous HST WFC3 spectrum (Benneke 
et al. 2019a) in the 1.1-1.7 um range, as shown in Fig. 6, 
there is notable difference in two of the data points at 
the blue end of the WFC3 spectrum. A comparison be- 
tween the new JWST spectrum and the HST spectrum is 
presented in the Appendix A. Furthermore, the presence 
of multiple CH, features across our NIRISS and NIR- 
Spec spectral range provide a very strong detection of 
CHA, as discussed in section 3.4. We note that our upper 
limit for H20 corresponds to the planet’s stratosphere 
at pressures below ~100 mbar. Water vapor may very 
well be abundant at deeper levels in the atmosphere, 
but condensation of H2O is expected in the upper tro- 
posphere of this temperate planet (Benneke et al. 2019a; 
Madhusudhan et al. 2023), resulting in a comparatively 
dry stratosphere, as on Earth. 

We also do not detect CO or HCN despite their strong 
spectral features expected in the 0.9-5.2 um range (Mad- 
husudhan et al. 2023). The 95% upper limits for both 
molecules are shown in Table 2, with maximum values of 
-3.00 for log(Xco) and -2.41 for log(Xucw). Given the 
low-temperature H5-rich atmosphere, the non-detection 
of CO is not necessarily surprising, as CH, is expected 
to be the dominant equilibrium constituent in deep H2 
atmospheres on cooler planets (Moses et al. 2013). How- 
ever, some CO is expected to be present from disequi- 
librium quenching in deep atmospheres or photochem- 
istry at high altitudes, becoming especially important in 
thinner atmospheres (Yu et al. 2021; Tsai et al. 2021; Hu 
et al. 2021; Madhusudhan et al. 2023). The high abun- 
dances of CO2 and CH,, along with the non-detection 
of NH; and CO, and a high CO2/CO ratio, are consis- 
tent with predictions for an ocean surface under a thin 
Ha-rich atmosphere (Hu et al. 2021; Madhusudhan et al. 
2023), as discussed further in section 4. 


3.3. Biosignature Molecules 


The retrievals provide notable constraints on two 
methyl-group terrestrial biomarkers, DMS and CH3Cl, 
predicted to be detectable in Hycean atmospheres, es- 
pecially for K2-18 b (Madhusudhan et al. 2021). We 


retrieve a DMS mixing ratio of log(Xpms) = —4.46*0 4 


in the no-offset case, —6.35*100 in the one-offset case 
and —6.87^1 57 in the two-offsets case. The weaker con- 
straints on the DMS abundance with increasing number 
of offsets between the detectors is due to the DMS spec- 
tral feature being broad and the continuum level span- 
ning multiple detectors, as discussed in section 3.4. The 
potential inference of DMS is of high importance as it is 
known to be a robust biomarker on the Earth and has 
been extensively advocated to be a promising biomarker 
for exoplanets (Seager et al. 2013; Seager et al. 2016; 
Catling et al. 2018; Madhusudhan et al. 2021); this is dis- 
cussed further in section 4. We also find a nominal peak 
in the posterior distribution of CH3Cl which is more sig- 
nificant than other non-detections such as those of HzO 
or NH3, as shown in Fig. 4. The retrieval in the no-offset 
case provides an abundance estimate of log(Xcn,ci) = 
—6.62*2 09 and a 95% upper limit of -2.50, which are 
comparable within ~1-dex with those from the other 
retrieval cases. 

We note that CH4, DMS and CH3CI all have strong 
spectral features due to the C-H bond in the 3-3.54m 
range in the NIRSpec G395H band and are therefore 
degenerate to some extent, as shown in Fig. 5. How- 
ever, thanks to the multiple strong CH, features in 
the NIRISS band, the degeneracies between the two 
stronger molecules, CH4 and DMS, are somewhat miti- 
gated, whereas CHCl is relatively unconstrained. 


3.4. Molecular Detection Significances 


In addition to the abundance constraints discussed 
above, we determine the detection significances of 
the key molecules using Bayesian model comparisons 
(Trotta 2008; Benneke & Seager 2013; Pinhas et al. 
2018). In the present context, we evaluate the detection 
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significance (DS) of a molecule as a Bayesian preference 
for a model fit to the data while including that molecule, 
relative to the same model with the molecule absent 
(Benneke & Seager 2013; Pinhas et al. 2018). Naturally, 
such an evidence comparison depends to some extent 
on the specific model considerations and combination of 
the data used, which are discussed in section 3.1. There- 
fore, we estimate detection significances for the promi- 
nent molecules for each of the three retrieval cases, as 
shown in Table 2. We note that the detection signifi- 
cances reported here have intrinsic statistical uncertain- 
ties of ~0.1o due to the uncertainty in Bayesian evidence 
obtained by the nested sampling algorithm for a given 
retrieval. 

Among the prominent CNO molecules, we find the 
strongest detection for CH, at 4.7-5.0c across all three 
cases. The consistently high DS value independent of 
the offset(s) considered underscores the robustness of 
the detection, which is due to multiple features of CH4 
being present across the 1-5 jum range of the observed 
spectrum as shown in Fig. 3. We also detect CO» ro- 
bustly at ~3ø significance across all the three cases. The 
strong detection of CO» is made possible by its promi- 
nent spectral feature around 4.3 jum and the full feature 
along with the spectral baseline around it being on the 
same detector (NRS2) in the NIRSpec band. As dis- 
cussed above, we do not find significant evidence for 
NHs, H20, CO or HCN. 

Among the biomarkers, we find some evidence for 
DMS depending on the retrieval case. The detection 
significance of DMS depends on the offsets considered. 
This is because, while DMS has a strong spectral fea- 
ture around 3.3 um in the NIRSpec NRS1 detector, the 
feature is broad and the spectral baseline falls on neigh- 
bouring detectors — the NIRISS at shorter wavelengths 
and NIRSpec NRS2 at longer wavelengths. Therefore, 
the spectral amplitude and hence the detection signif- 
icance and abundance estimate relies strongly on the 
relative offsets between the detectors; the detection sig- 
nificance lowers for each additional offset in the retrieval. 
We infer DMS at 2.40 confidence for the no-offset case, 
but at only ~1o for the 1-offset case and no significant 
evidence for the 2-offset case. Nevertheless, as shown in 
Fig. 4, in all three cases the retrieved posterior distribu- 
tions for DMS show notable peaks within 1-dex of each 
other, except that for the cases with offsets the distribu- 
tions contain long low-abundance tails owing to the de- 
generacy with the spectral baseline as discussed above. 
The posteriors are also notably different from the non- 
detections for other prominent molecules such as H20, 
NH; or HCN, and indicate a non-negligible chance for 
DMS being present in the atmosphere. Upcoming ob- 


servations will be able to further constrain the presence 
of DMS, as discussed below and in section 4. 

There could also be potential contributions from 
CH3Cl to the spectrum, albeit without any appreciable 
detection significance, as evident from Fig. 4. CH3Cl, 
being a methylated molecule like DMS, has some over- 
lapping features with DMS as shown in Fig. 5. There- 
fore, its significance increases marginally in the absence 
of DMS. While we do not detect CH3Cl on its own in 
any of the retrieval cases, the combination of DMS and 
CH3CI has a slightly higher detection significance of 2.70 
than DMS alone (2.40) in the no-offset case. We find no 
significant evidence for any other biomarkers considered 
in the retrievals. 

Overall, we find CH4 and CO» to be our most confi- 
dent detections, followed by DMS, with the abundance 
estimates reported above. While our results provide im- 
portant first insights into the chemical composition of 
K2-18 b, upcoming observations will be able to verify 
our present findings. These include observations of the 
transmission spectrum of K2-18 b with JWST MIRI be- 
tween ~5-10 um (JWST Program GO 2722, PI: N. Mad- 
husudhan) and more observations with JWST NIRSpec 
G395H and G235H (JWST Program GO 2372, PI: R. 
Hu). 


3.5. Clouds/Hazes and Photospheric Temperature 


'The observed transmission spectrum provides nominal 
constraints on the presence of clouds/hazes in the atmo- 
sphere. The constraints on the cloud/haze parameters 
are shown in Table 3. The constraints on the cloud-top 
pressure for the gray clouds are relatively weak, mostly 
lying below the observable photosphere (e.g. cloud top 
pressures 2100 mbar). Even though the scattering slope 
(y) is not well constrained, the enhancement factors (a) 
are generally higher than that expected for Hg Rayleigh 
scattering (a—1), albeit still consistent with the latter 
at the 3e uncertainties for the no-offset case. The haze 
coverage fraction at the day-night terminator region is 
constrained to ~0.6, albeit with large uncertainties of 
~0.2. Based on Bayesian model comparisons as dis- 
cussed above, we find that a model with clouds/hazes is 
preferred over a model without clouds/hazes by 2.8-3.20 
across the three retrieval cases considered. However, 
more observations in the optical to near-infrared wave- 
lengths would be needed to provide stronger constraints 
on the cloud/haze properties of K2-18 b, as discussed in 
Appendix B. We find that the abundance constraints on 
the molecules are consistent within the 1-o uncertainties 
between the retrievals with and without clouds/hazes. 

We additionally consider retrievals in which the para- 
metric clouds/hazes of the canonical model are replaced 
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with Mie scattering hazes, as described by Pinhas & 
Madhusudhan (2017) and Constantinou et al. (2023). 
We specifically include two forms of organic haze, using 
optical constants presented by Khare et al. (1984) and 
He et al. (2023). We find no evidence for this model, 
and a Bayesian model comparison shows a preference in 
favour of the parametric clouds/hazes considered in our 
canonical model. While the haze properties are uncon- 
strained in this retrieval, the abundance constraints on 
the gaseous species remain consistent with those from 
our canonical retrieval cases. 

'The observations provide nominal constraints on the 
temperature in the planetary photosphere. We retrieve 
the temperature at 10 mbar to range between 235728 
and 257*17' between the three cases, as shown in Ta- 
ble 3. We note that the retrieved temperature struc- 
ture is typically less well constrained using transmis- 
sion spectroscopy, compared to emission spectroscopy 
(Madhusudhan et al. 2016). Nevertheless, the retrieved 
temperature range and the non-detection of H20 allow 
for the possibility of H20 clouds in the deeper atmo- 
sphere. Considering the pressures probed by the spec- 
tral features across our observed range, we find that the 
photosphere, ie. the 7—1 surface, lies between pres- 
sures of ~0.1-100 mbar. The 10 mbar temperature esti- 
mates for the three retrieval cases are shown in Table 3. 
While the H2-rich atmosphere can result in a significant 
greenhouse effect and warm the ocean surface, clouds 
and/or hazes play a crucial role in cooling the atmo- 
sphere and decreasing the temperature gradient (Mad- 
husudhan et al. 2020, 2021, 2023; Piette & Madhusud- 
han 2020). The possible presence of clouds can allow 
more temperate conditions at the ocean surface com- 
pared to those predicted by cloud-free models (cf. Innes 
et al. 2023). 


3.6. Stellar Heterogeneities 


We also consider retrievals including the effects of 
unocculted stellar heterogeneities on the transmission 
spectrum, using our AURA retrieval framework (Pin- 
has et al. 2018). We do not find significant evidence for 
such effects in any of the three retrieval cases. The spot- 
covering fraction in our retrieval is consistent with zero 
at the 2e uncertainties, and a model with stellar het- 
erogeneity is not favored over a model without it across 
the three cases. We also find that the abundance con- 
straints on the molecules are not significantly affected 
by the consideration of stellar heterogeneities in the re- 
trieval. The non-detection of H20 is further evidence 
against the effects of unocculted stellar heterogeneities 
on the present transmission spectrum, given that H20 
in cool unocculted starspots may contaminate the spec- 


trum, as recently reported (Barclay et al. 2021; Moran 
et al. 2023). 


4. SUMMARY AND DISCUSSION 


We report a transmission spectrum of the candidate 
Hycean exoplanet K2-18 b observed with JWST. The 
spectrum observed with the JWST NIRISS and NIR- 
Spec instruments spans the 0.9-5.2 um range containing 
strong absorption features of prominent CNO molecules 
and biomarkers predicted for Hycean worlds. We re- 
port strong detections of CH, and CO»; in a Hə-rich 
atmosphere at 5g and 30 confidence, respectively, with 
high volume mixing ratios (~1%) for both. However, 
we do not detect H20, NH3, CO or HCN, while ob- 
taining upper-limits on their abundances that are con- 
sistent with chemical expectations for an ocean under a 
cold and thin H5-rich atmosphere (Hu et al. 2021; Mad- 
husudhan et al. 2023). We also find potential evidence 
for dimethyl sulfide (DMS), which has been predicted 
as a robust biomarker in both terrestrial and Hycean 
worlds. These findings support the Hycean nature of 
K2-18b and the potential for biological activity on the 
planet. 

The observed mass, radius and equilibrium tempera- 
ture of K2-18b have been known to be consistent with 
a degenerate set of internal structures (Madhusudhan 
et al. 2020, 2023). These include (a) a Hycean world 
with a thin Hə-rich atmosphere over a water-rich in- 
terior, (b) a mini-Neptune with a deep Ho-rich atmo- 
sphere, or (c) a predominantly-rocky super-Earth in- 
terior with a deep Ho-rich atmosphere. Our retrieved 
chemical composition of the atmosphere of K2-18b helps 
distinguish between these scenarios. In what follows, we 
discuss the implications and possible explanations of our 
findings and future directions. 


4.1. A Potential Hycean World 


K2-18 b has been originally predicted to be the 
archetype of a Hycean world (Madhusudhan et al. 2021), 
one with habitable oceans underneath a Hə-rich atmo- 
sphere. The presently derived chemical composition of 
the atmosphere is in agreement with previous theoreti- 
cal predictions for the presence of an ocean under a shal- 
low Ha-rich atmosphere (Hu et al. 2021; Madhusudhan 
et al. 2023). In particular, Hu et al. (2021) predicted 
the abundance of CO» to range between 4x10~* and 
107} and that of CH, between 1.5-5.3 x1072, which are 
in agreement with our retrieved abundances. They also 
predicted relatively lower abundances of CO, NH3 and 
stratospheric H20, which are consistent with our non- 
detections. A low retrieved H20 gas-phase mixing ratio 
at pressures less than ~100 mbar is consistent with con- 
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Figure 5. Spectral contributions of key molecular species in the 1-5 um range. The different curves show individual contributions 
from different molecules to a nominal model transmission spectrum of K2-18 b shown in blue and denoted as Combined. 
The model assumes a mixing ratio of 107? for CH4 and CO», 1074 for H2O, 10^? for all the other species, consistent with 
our retrieval estimates discussed in section 3, and an isothermal temperature profile of 250 K. Each curve corresponds to a 
transmission spectrum with opacity contributions from a single molecule at a time, in addition to H2-Hz and H»-He collision- 
induced absorption. The spectral ranges of our JWST NIRISS and NIRSPec observations are also indicated; the NIRSpec range 
spans two detectors (NRS1 and NRS2) with a gap between them at 3.72-3.82 um. 


densation due to a tropospheric cold trap (Madhusud- 
han et al. 2023), as in the Earth's stratosphere, and with 
the retrieved thermal structure in this work and previ- 
ous studies (e.g. Benneke et al. 2019a; Piette & Mad- 
husudhan 2020; Madhusudhan et al. 2023). That is, 
H20 could be abundant below its condensation region 
in the atmosphere, but the transit observations of the 
terminator region do not probe deep enough to detect 
it. 

The main argument against K2-18 b being a poten- 
tial Hycean world is based on climate considerations — 
a greenhouse effect in a thick H5-rich atmosphere that 
is cloud/haze-free would result in temperatures being 
sufficiently elevated at pressures greater than ~10 bar, 
such that a liquid ocean would instead be converted 
to a steam-dominated atmosphere that ultimately goes 
supercritical at depth for irradiation levels relevant to 
K2-18 b (Scheucher et al. 2020; Piette & Madhusud- 
han 2020; Innes et al. 2023; Pierrehumbert 2023). Any 
global ocean surface must then reside at pressures less 
than ^10 bar. However, too shallow a Hz atmosphere 
could be subject to escape over time (e.g., Kubyshkina 
et al. 2018a,b; Hu et al. 2023), so there is a limited pa- 
rameter range over which a habitable ocean could exist 
on K2-18 b without significant clouds/hazes. As men- 
tioned previously, high-albedo tropospheric water clouds 
or scattering hazes can help alleviate the steam and su- 
percritical water problem by reducing the stellar energy 


absorbed by the planet (Madhusudhan et al. 2021; Piette 
& Madhusudhan 2020). 

Our retrieved atmospheric temperatures and upper 
limits for H20 are consistent with the possibility that 
H20 is condensing into clouds below the photosphere 
in K2-18 b, indicating a cold upper troposphere (see 
also Benneke et al. 2019a; Madhusudhan et al. 2023). 
Furthermore, as discussed in section 3.5, our retrievals 
show some evidence for scattering due to hazes at the 
day-night terminator region of the atmosphere; though, 
more optical observations are required to robustly con- 
firm the same. If such clouds or hazes enshroud the 
planet or, in particular, are present on the dayside, they 
could provide the required albedo to sustain a habitable 
ocean in K2-18 b. Given the possible greenhouse effects 
of clouds themselves, a relatively shallow H5 atmosphere 
might still be required to maintain low-enough temper- 
atures for a liquid ocean. 


4.2. Is a Deep Atmosphere a Possibility? 


Our derived atmospheric composition in combination 
with K2-18 b's irradiation level is seemingly inconsistent 
with a deep-atmosphere, mini-Neptune scenario. In that 
scenario any photochemically produced carbon and ni- 
trogen species would be recycled in the deep atmosphere 
back to their thermodynamically stable forms, CH4 and 
NHs, with transport then returning these molecules 
back to the observable upper regions of the atmosphere 
(Yu et al. 2021; Tsai et al. 2021; Madhusudhan et al. 
2023). Although CO» can have a substantial mixing 
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ratio in a Ha5-rich atmosphere with a high metallicity 
and/or low C/O ratio (e.g., Moses et al. 2013), we find 
that the simultaneous presence of ~1% COs and CH, in 
a thick H5-rich atmosphere requires a moderately high 
C/H metallicity, a very low C/O ratio, (70.02) and ef- 
ficient vertical quenching. For example, C/H — 30x 
solar, O/H = 690x solar, and K,, Z 10" cm? s^! com- 
bined with a temperature profile with Tine = 60 K pro- 
vides CH4 and CO» mixing ratios of the right magni- 
tude. Aside from the question of how such an oxygen- 
enriched atmosphere would originate from a formation 
and evolution standpoint, this scenario would cause H20 
to dominate over Hə throughout much of the atmosphere 
below a few hundred mbar (increasing its mean molec- 
ular weight), and NH3 and CO would also be present in 
non-negligible quantities, none of which are consistent 
with the observations. 

Nor can our observations be explained by a shal- 
low, Ha-rich atmosphere overlying a solid surface at the 
few-bar pressure level, despite the abundances of CH4 
and CO» being seemingly consistent with such models 
(Yu et al. 2021; Tsai et al. 2021; Madhusudhan et al. 
2023). The planet's bulk density precludes a thin Hə 
atmosphere overlying an extensive silicate mantle (Mad- 
husudhan et al. 2020). Even a pure silicate interior 
would require a thick (Z10? bar) Ha-rich envelope to 
explain the mass and radius. Although interaction of a 
thick Hz atmosphere with a deep silicate mantle might 
explain the presence of some atmospheric CO» (e.g., 
Kite et al. 2020; Schlichting & Young 2022; Tian & 
Heng 2023), such models do not predict a resulting 196 
CO, mixing ratio from high-pressure interaction with 
silicates at depth in a thick Hz atmosphere, nor can 
that scenario explain the apparent lack of recycling of 
Nə back to NH3. Overall, the bulk density of the planet 
combined with our derived chemical composition for the 
atmosphere present strong evidence in support of K2- 
18 b as a Hycean world rather than a rocky or volatile- 
rich planet with a deep Hz atmosphere, or a rocky planet 
with a thin H9 atmosphere. 

One caveat to the above discussion is that current pho- 
tochemical models for mini-Neptune conditions assume 
ideal-gas behavior and do not consider how a primor- 
dial H5 atmosphere with its expected reduced species 
might interact chemically with a supercritical water 
layer and/or silicate magma at depth in the atmosphere. 
Moreover, our understanding of the bulk composition 
and chemistry of super-Earth and mini-Neptune atmo- 
spheres is rudimentary at this stage, and the depleted 
NH; and CO could potentially have some chemical ex- 
planation that has not been considered up to this point. 
We also note that the retrieved abundances represent 


average chemical abundances in the observable photo- 
sphere (~0.1-100 mbar) at the day-night terminator re- 
gion. While the molecules with robust detections, CH4 
and CO», have been predicted to be relatively uniform 
in this pressure range for K2-18 b in several cases (Hu 
et al. 2021; Yu et al. 2021; Madhusudhan et al. 2023) 
more precise observations in the future may be able to 
constrain non-uniform chemical abundances in the at- 
mosphere. 


4.3. Possible Evidence of Life 


Our potential evidence for DMS in K2-18 b moti- 
vates consideration of possible biological activity on the 
planet. While the present evidence is not as strong as 
that for CH4 or CO2, upcoming JWST observations of 
K2-18 b will be able to robustly constrain the presence 
and abundance of DMS, as discussed in section 4.5 and 
earlier work (Madhusudhan et al. 2021). Here we dis- 
cuss the plausibility of our DMS abundance constraints 
from a potential biosphere on K2-18 b in order to inform 
future observations and retrieval studies. 

On Earth, DMS is a byproduct of living organisms, 
with the bulk of the non-anthropogenic DMS in the ter- 
restrial atmosphere being emitted from phytoplankton 
in marine environments (Charlson et al. 1987; Barnes 
et al. 2006). Both DMS and CH3Cl are thought to be 
terrestrial biosignatures with no known false positives 
(Catling et al. 2018). On the Earth these molecules are 
produced exclusively by life in relatively small quanti- 
ties compared to more abundant byproducts of life, such 
as O2, CH4 and N20; the latter are, therefore, more 
favoured as biosignatures for Earth-like planets. How- 
ever, it has been suggested that in Hg-rich environments 
with large biomass, molecules such as DMS and CH3Cl 
could be abundant and observable for habitable super- 
Earths (Seager et al. 2013) and Hycean worlds (Mad- 
husudhan et al. 2021). 

While we infer DMS with marginal confidence, our 
retrieved DMS abundance spans a relatively wide 
range across the cases considered, from log(Xpms) = 
—4.46*0 1. at the higher end for the no-offset case to 
log(Xpms) = —6.87*157 for the two-offset case. On 
Earth, the typical mixing ratios of DMS, found near 
the ocean surface, are a few hundred parts per trillion 
(Hopkins et al. 2023). DMS is rapidly depleted at higher 
altitudes, with a few-day lifetime, due to photochemical 
reactions with OH and other radicals, ultimately lead- 
ing to the production of more oxidized sulfur molecules, 
such as SO». The upper end of our retrieved abundance 
for DMS is significantly higher than that on Earth, and 
we do not detect other sulfur-bearing species; however, 
the lower end is more plausible (Seager et al. 2013). We 
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note that the infrared absorption cross-sections of DMS 
are currently limited (Kochanov et al. 2019; Sharpe et al. 
2004; Gordon et al. 2017). New cross-section data may 
revise our mixing ratio estimates, while future more in- 
tensive observations could be sensitive to other sulfur 
species. 

'The DMS abundance is also strongly dependent on the 
chromospheric activity of the host star. A quiescent M 
dwarf host star with a lower ultraviolet flux compared to 
a sun-like star of the same bolometric flux could enable 
DMS to survive longer and be more abundant in the 
planetary atmosphere (Domagal-Goldman et al. 2011). 
Previous studies have predicted DMS mixing ratios of 
7.10-7—10-5 for Earth-like planets and super-Earths or- 
biting low-activity M dwarfs for plausible biomass es- 
timates (Domagal-Goldman et al. 2011; Seager et al. 
2013). However, K2-18 b has been described as a mod- 
erately active M3 dwarf (Benneke et al. 2019b), and it 
may not be quiescent in the extreme ultraviolet (dos 
Santos et al. 2020). A high Lyman alpha flux could 
lead to DMS depletion due to interaction with atomic 
O produced from CO» photolysis (Seager et al. 2013). 
Reaction of DMS with atomic H is significantly slower 
than with O (cf. Atkinson et al. 2004; Zhang et al. 2005), 
and the OH abundance would be limited by H20 con- 
densation in K2-18 b. 

If the DMS abundance on K2-18 b is indeed confirmed 
by future observations to be greater than ~10~°, that re- 
sult could require very high biological production rates 
in the ocean and/or new theoretical developments in 
our understanding of DMS chemistry (including poten- 
tial abiotic chemistry) in planets such as K2-18 b. We 
also note that while our retrievals included a wide array 
of molecules with strong spectral signatures in the ob- 
served wavelength range, future theoretical studies and 
retrievals could consider an even more expanded set, 
particularly as accurate cross-section data become avail- 
able. In particular, other hydrocarbon molecules have 
similar C-H bands in the 3-3.5 um range as DMS. There- 
fore, besides our canonical model, we have explored re- 
trievals with other hydrocarbons such as C2H2, C2He, 
CH3OH, HC3N, and hazes with optical properties from 
He et al. (2023), which have been predicted to be rele- 
vant for Hycean atmospheres (e.g. Tsai et al. 2021; Mad- 
husudhan et al. 2023), but we found no definitive evi- 
dence for their presence. 

Besides DMS, the presence of life could potentially 
also contribute to the strong chemical disequilibrium in- 
dicated by the retrieved atmospheric composition of K2- 
18 b. For instance, methanogenic bacteria in Earth's 
oceans are known to be a significant contributor to the 
atmospheric CH4 budget. It is possible that similar bi- 


otic sources may to some extent contribute to the ob- 
served CH4 abundance in K2-18 b, if indeed life exists 
on the planet. 

Overall, our findings demonstrate the feasibility of de- 
tecting a biosignature molecule in the atmosphere of a 
habitable-zone sub-Neptune with JWST. This also pro- 
vides a valuable case study for a framework for biosigna- 
ture assessment in exoplanets (e.g. Catling et al. 2018; 
Meadows et al. 2022). The potential inference of DMS 
in K2-18 b provides a pathway towards the possible de- 
tection of life on an exoplanet with JWST and other 
current and upcoming large observational facilities. The 
next steps would involve both (a) more theoretical in- 
vestigations to understand the possible atmospheric and 
interior processes at play, and (b) more observations to 
verify the present findings and potentially discover other 
chemical species. 


4.4. Resolving the Missing Methane Problem 


Our strong detection of CH4 at 5e resolves one of 
the longest standing conundrums in exoplanet science 
— “The Missing Methane Problem" (Stevenson et al. 
2010; Madhusudhan & Seager 2011). Low-temperature 
molecules such as CH4 and NH3 are common in the so- 
lar system and are seen in the atmospheres of the giant 
planets (Karkoschka 1998; Encrenaz et al. 2005; Atreya 
et al. 2018). These molecules are expected to be promi- 
nent carriers of carbon and nitrogen in Ho-rich atmo- 
spheres at temperatures below ~600 K, with H20 being 
the dominant oxygen carrier (Burrows & Sharp 1999; 
Lodders & Fegley 2002). However, no robust detection 
of CH4 or NHs has been made in any exoplanetary at- 
mosphere with temperatures below ~800 K, despite at- 
mospheric observations made for several such exoplanets 
with HST and Spitzer at wavelengths sensitive to these 
molecules, e.g. GJ 436 b (Stevenson et al. 2010; Knutson 
et al. 2014), GJ 3470 b (Benneke et al. 2019b), K2-18 b 
(Benneke et al. 20192), all in the ~300-800 K range; (cf. 
Blain et al. 2021). 

Atmospheres of temperate sub-Neptunes are expected 
to exhibit distinct non-equilibrium chemistry, just as in 
solar system planets. Several processes can cause chem- 
ical disequilibrium, including photochemistry, vertical 
mixing and volcanic outgassing (e.g. Yung 1999). How- 
ever, even strong non-equilibrium chemical mechanisms 
have difficulty explaining the missing CH, in temper- 
ate exoplanetary atmospheres (Line et al. 2011; Moses 
et al. 2013). This missing methane problem has, there- 
fore, remained one of the central puzzles in the area of 
exoplanetary atmospheres in the pre-JWST era. The 
present detection of CH4 in K2-18 b, therefore, demon- 
strates the detectability of CH4 and potentially other 
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hydrocarbons with JWST and opens a new era of atmo- 
spheric characterisation of temperate exoplanets in gen- 
eral. The detection suggests that other sub-Neptunes 
and giant exoplanets with Ho-rich atmospheres at sim- 
ilar temperatures may also be conducive for detecting 
CH,4. Such planets therefore represent important tar- 
gets for homogeneous studies of carbon chemistry in ex- 
oplanetary atmospheres, enabling comparative studies 
with solar system giant planets. 


4.5. Future Directions 


Our results demonstrate the potential of candidate 
Hycean worlds as optimal targets in the search for life 
on exoplanets. These findings motivate further obser- 
vations and theoretical work to characterise in detail 
the atmospheric and potential surface conditions of K2- 
18 b and other candidate Hycean worlds (Madhusudhan 
et al. 2021). Several upcoming JWST observations of 
K2-18 b will be able to verify the present findings, in par- 
ticular more observations with NIRSpec G395H (JWST 
GO 2372) and MIRI LRS (5-10 um) (JWST GO 2722). 
While the former program can confirm the present find- 
ings with higher precision, the latter can specifically con- 
firm the presence of DMS which is expected to have a 
larger feature around 7 um (Seager et al. 2013). Such 
observations are also motivated for a number of other 
promising candidate Hycean worlds orbiting nearby M 
dwarfs that are even more favorable to observations than 
K2-18 b (Madhusudhan et al. 2021). 

Overall, the present results pave the way to a new 
era of atmospheric characterisation of habitable planets 
and biosignature detection with JWST. The observa- 
tions also motivate a wide range of theoretical studies to 
understand in detail the physical, chemical, and biologi- 
cal conditions on Hycean worlds. Our findings present a 
first step towards the spectroscopic identification of life 
beyond the solar system and the assessment of our place 
in the Universe. 
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APPENDIX 


A. COMPARISON WITH HST WFC3 OBSERVATIONS 


As discussed above, K2-18 b has been previously observed in the 1.1-1.7um range with HST WFC3 (Tsiaras et al. 
2019; Benneke et al. 2019a), with these observations leading to inferences of H20 being present in the planet's ter- 
minator atmosphere (Tsiaras et al. 2019; Benneke et al. 2019a; Madhusudhan et al. 2020) with a high upper-limit on 
CH, (Madhusudhan et al. 2020). The spectrum has also been suggested to be explained by CH4 or a combination of 
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Figure 6. HST and JWST observations of K2-18 b. The black points show prior observations of K2-18 b obtained with HST 
WFC3 in the 1.1-1.7 um range, with data from 8 transits, presented by Benneke et al. (2019a) using data from HST GO 13665 
and GO 14682 programs (PI: B. Benneke). The orange and dark red points show our JWST NIRISS and NIRSpec observations 
from one transit each. The data is binned for visual clarity to R = 25 and R = 55, respectively, as shown in Fig. 3. The dark 
blue line denotes the median retrieved spectrum (one-offset case), while medium and lighter blue regions denote the lo and 26 
contours, respectively. The yellow points correspond to the median spectrum binned to match the JWST observations. Our 
JWST NIRISS spectrum is in agreement with the HST WFC3 spectrum for most of the common wavelength range except for 
two data points towards the blue end of the WFC3 band. 


CH, and H20, instead of H20 alone, due to the degeneracy between the two molecules in the WFC3 band (Blain et al. 
2021). Fig. 6 overlays these prior observations presented by Benneke et al. (2019a) with our new JWST spectrum. It 
can be seen that the new NIRISS data presented in this work is in good agreement with the prior HST WFC3 data in 
general, with the notable exception of two data points in the blue end of the WFC3 band, which show a 2-30 deviation 
from the NIRISS observations and corresponding spectral fit. The two deviant WFC3 points are inconsistent with a 
CH, absorption peak in the retrieved spectral fit. As such, it is possible that prior inferences of HzO over CH4 were 
driven by these two points, as H20 does not display a prominent absorption feature at those wavelengths 


B. COMPARISON WITH RETRIEVALS ASSUMING NO CLOUDS 


As discussed in section 3.5 our retrievals with the canonical model provide nominal constraints on the properties 
of possible clouds/hazes. Here we show a retrieved spectral fit to our JWST transmission spectrum using a model 
with no clouds/hazes (Fig. 7). This is applied to the 1-offset case discussed in section 3, with all other factors being 
the same as for the canonical retrieval with clouds/hazes included. While the Bayesian evidence is higher for the 
canonical model, as discussed in section 3.5, the no clouds/hazes model provides a reasonable fit (Fig. 7) to most 
of the observed spectrum. It differs from the canonical model fit (Fig. 3) only towards the blue end of the NIRISS 
spectrum. The retrieved abundances in the present case are also consistent with the canonical case within the lo 
uncertainties. Therefore, further observations in the optical are needed to more robustly constrain the presence and 
properties of clouds/hazes in the atmosphere of K2-18 b. 


C. ATMOSPHERIC RETRIEVAL BAYESIAN PRIORS 


Table 4 shows the Bayesian prior probability distributions used in the retrievals presented in this work. All but the 
last six parameters correspond to the canonical model described in section 3.1. dNirspec denotes the linear offset applied 
to the NIRSpec G395H data in the 1-offset case in parts-per-million, while ônrsı and Óxngo denote the separate offsets 
applied to observations from the NIRSpec NRS1 and NRS2 detectors, respectively, in parts-per-million. Tphot; Thet 
and fret are the photospheric temperature, starspot temperature and coverage fraction, respectively, and are used 
in the unocculted stellar heterogeneities modelling discussed in section 3.6. U(X,Y) denotes a uniform probability 
distribution between X and Y, while A (u,a?) denotes a normal distribution of mean u and variance o?. For any 
additional chemical species considered, such as the hydrocarbons described in section 4.3, we use the same uniform 
mixing ratio priors as those shown in table 4. 
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Figure 7. Retrieved spectral fit of the JWST transmission spectrum of K2-18b using a model without clouds/hazes in the 
one-offset case. Details of the figure are the same as for Fig. 3 except that the retrieval does not include clouds/hazes in the 
model. The retrieval produces a comparable fit to the cloudy case (Fig. 3) across most of the wavelength range, differing mostly 
towards the bluer end of the NIRISS band, as discussed in Appendix B. The NIRSpec spectrum is vertically offset by —30 ppm, 
corresponding to the median retrieved offset. 
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Description 


Parameter Bayesian Prior 
log(Xu,0) U(-12, -0.3) 
log(Xcu, ) U(-12, -0.3) 
log( Xwnuz ) U(-12, -0.3) 
log(Xucn) U(-12, -0.3) 
log(Xco) U (-12, -0.3) 
log(Xco;) U(-12, -0.3) 
log(Xpms) U(-12, -0.3) 
log(Xcs,) U(-12, -0.3) 
log(Xcu3c1) U (-12, -0.3) 
log(Xocs) U(-12, -0.3) 
log(Xn,0) U (-12, -0.3) 
To/K U(0, 500) 
o;/K^3 (0.02, 2.00) 
a /K~2 U(0.02, 2.00) 
log(P; /bar) U(-6, 0) 
log(P2/bar) U(-6, 0) 
log(P3/bar) U(-2, 0) 
log (Pret /bar) U(-6, 0) 
log(a) U (-4, 10) 
2 u (-20, 2) 
log( P. /bar) U(-6, 1) 

$ u(0, 1) 
ONirspec/ppm | U(-100, 100) 
ÔNRS1/ppm U(-100, 100) 
ÓNRs2/ppm U(-100, 100) 
Tonoc/ K N (3457, 100?) 
Thet/ K U (1729, 4148) 
net U(0, 0.5) 


Mixing ratio of H20 
Mixing ratio of CH4 
Mixing ratio of NH3 
Mixing ratio of HCN 
Mixing ratio of CO 
Mixing ratio of CO2 
Mixing ratio of DMS 
Mixing ratio CS2 
Mixing ratio of CH3Cl 
Mixing ratio of OCS 
Mixing ratio N20 


Temperature at lubar 

P-T profile curvature 

P-T profile curvature 

P-T profile region limit 

P-T profile region limit 

P-T profile region limit 
Reference pressure at Rp 
Rayleigh enhancement factor 
Scattering slope 

Cloud top pressure 
Cloud/haze coverage fraction 
NIRSpec dataset offset 
NIRSpec NRS1 dataset offset 
NIRSpec NRS2 dataset offset 
Stellar photosphere temperature 
Starspot temperature 


Starspot coverage fraction 


Table 4. Model parameters and corresponding Bayesian prior probability distributions used in the atmospheric retrievals 
performed in this work. All but the last six parameters are used in the canonical atmospheric model described in section 3.1. 
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